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^r^ ■ density of a quark-gluon plasma are presented. These approximations incorporate 

^D \ the essential physics of the hard thermal loops, involve only ultraviolet-finite quan- 

tities, and are free from overcounting ambiguities. While being nonperturbative in 
Q^ ■ the strong coupling constant (7, agreement with ordinary perturbation theory is 

achieved up to and including order c^ . It is shown how the pressure can be re- 
CIhI constructed from the entropy and the baryon density taking into account the scale 

Q.' anomaly. The results obtained are in good agreement with available lattice data 

D , down to temperatures of about twice the critical temperature. 
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Self-consistent approximations allowing the calculation of the entropy and the baryon 



1 Introduction 



While asymptotic freedom suggests that, at very high temperature and/or suf- 
ficiently large baryon density, QCD should behave as a weakly coupled quark- 
gluon plasma, lattice results reveal that this ideal gas limit is approached 
only rather slowly. Perturbation theory [1,2] fails to reproduce this behaviour 
correctly, and resummations beyond the conventional ring-resummation are 
needed to properly account for the thermodynamics. In QCD, the situation 
may not be hopelessly complicated however; in fact, the available lattice data 
are quite well reproduced by simple models based on non-interacting "quasi- 
particles" with temperature dependent masses [3,4], which suggests that it 
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could be possible to describe the thermodynamics of the quark-gluon plasma 
in terms of its elementary excitations. 

However, phenomenological fits using massive quasiparticles represent only 
a convenient parametrisation of the thermodynamics and further theoretical 
effort is needed to arrive at a complete understanding in terms of the elemen- 
tary degrees of freedom and their interactions. What the quasiparticle picture 
suggests is that much of the physics at high temperature can presumably be 
described by propagator renormalisation or, in other words, a resummed per- 
turbation theory in terms of dressed propagators. The importance of retaining 
their effects nonperturbatively has been demonstrated recently in scalar toy 
models [5,6]. The general framework for carrying out consistently such renor- 
malisations of the thermodynamical functions has been well developed in the 
past in the context of the non-relativistic many-body problem [7-9] . 

A convenient quantity to consider when doing propagator renormalisation is 
the entropy, rather than the pressure [9,10]. Entropy is better suited, indeed, 
because it has a more direct interpretation in terms of quasiparticles; in par- 
ticular, it is expressible entirely in terms of dressed propagators, at least in 
lowest orders in a skeleton expansion. Building upon such an expansion, we 
are able to construct approximations which incorporate, in a way that is free 
of overcounting ambiguities, all the important physics of the dominant contri- 
butions to the self-energies at high temperature, the so-called "hard thermal 
loops" (HTL) [11,12]. (Corresponding vertex corrections, a priori needed to 
maintain gauge symmetry, turn out to be unimportant at the order of in- 
terest.) In contrast to the direct HTL resummation of the pressure recently 
proposed in Ref. [13], the resummation that we perform involves only UV- 
finite expressions and when compared to perturbation theory it agrees up to, 
and including, order g^. 

In Ref. [14], this strategy has been successfully applied to estimate the entropy 
S of the purely Yang-Mills plasma. In this Letter, we consider two important 
extensions. The first one is the inclusion of quarks which allows us to also treat 
plasmas with a finite chemical potential /i. The relevant formalism, presented 
in Sees. 2 and 3 below, leads to self-consistent approximations for both the 
entropy S and the baryon density Af. The second extension concerns the 
computation of thermodynamical quantities other than the entropy or the 
density. In Sec. 4, we show how to reconstruct the pressure P by a numerical 
integration of S and A/"; because of the trace anomaly, this requires some 
nonperturbative input which we take from lattice simulations at /i = 0. For 
all temperatures T ^ 2 . . . 3Tc, our results are in remarkable agreement with 
the corresponding lattice results whenever the latter are available. We also 
present results for the baryon density at T = and large //. Further possible 
extensions of our method are outlined in Sec. 5. 



2 Self-consistent approximation to the entropy and quark density 



The thermodynamic potential Q = —PV of a field theory involving bosons 
and fermions, expressed as a functional of full propagators {D for bosons, S 
for fermions) has the form [7] 



n[D, S] = |TTr logD-i - iTTr HD - TTr logS"^ + TTr ES + T$[D, S] 

= T^D,S]+ tr f -^^n{iu)lm [\ogD'\u,k) -n{to,k)D{u,k) 



d^k 



+2 try -r-_/(cu)Im [log5-i(cc.,A;)-S(cu,A;)5(cu,A;) 

where $[-D, S*] is the sum of 2-particle- irreducible "skeleton" diagrams, n{uj) = 
{e^^-l)-\ f{uj) = (e^('^-'^)+l)-i, and f3 = 1/T. Here "tr" refers to all dis- 
crete labels, including colour and fiavour when applicable. 

The self-energies II = D^^ — Dq^ and T, = S^^ — Sq'^, where Dq and 5*0 are 
bare propagators, are themselves functionals of the full propagators and are 
determined by 

S<!>[D,S]/dD = lU, S<^[D,S]/SS=J:, (2) 



which goes hand in hand with the all-important stationarity property 

6Q[D,S]/6D = = 6Q[D,S]/6S. (3) 

A "self-consistent" ( "^-derivable" ) [8] approximation is one that preserves this 
stationarity property by selecting a subset of skeleton-diagrams from $ and 
determining the self-energies from (2). In particular, a two- loop approximation 
to $[D,S'] (corresponding to a dressed one-loop approximation for the self- 
energies) is obtained by discarding all skeleton diagrams of loop-order 3 and 
higher. 

This two-loop approximation for Q has a remarkable consequence for the first 
derivatives of the thermodynamic potential, the entropy and the fermion den- 
sities: 
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Because of the stationarity property (3), one can ignore the T and /i depen- 
dences implicit in the spectral densities of the full propagators, and differen- 
tiate exclusively the statistical distribution functions n and / in (1). Now the 



derivative of the two-loop functional T^[D, S] at fixed spectral densities of the 
propagators D and S turns out to just cancel that part of the terms Im (JID) 
and Im (SS*) in (1) which involves RelllmD and Re Sims', respectively. 
In a self- consistent two-loop approximation one therefore has the remarkably 
simple formulae 



S = -ti J -r^^-^[lm\ogD-\u;,k)-lmU{u;,k)ReD{u;,k) 

-2tT J -—--!^[lm\ogS-\u,k)-lmE{u,k)ReS{u,k)],{5) 
A/'=-2tr / , ,, ■'^ ' \m\ogS-\uj,k) - \mT.{uj,k)ReS{uj,k) .(6) 
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This has been noted first for the entropy in a non-relativistic context by Riedel 
[10] and more recently in QED by Vanderheyden and Baym [15], but this 
important simplification holds much more generally and is equally valid for 
the fermion density [16]. 

In a nonabelian gauge theory, the above relations have to be in general aug- 
mented by ghost contributions, which can however be avoided in those gauges 
where ghosts do not propagate, such as Coulomb gauge, in which we shall work 
in what follows. At any rate, in gauge theories, abelian as well as nonabelian, 
self-consistency does not guarantee gauge invariance, for only propagators have 
been dressed, and no vertices^. 

However, in what follows we shall consider approximately self- consistent ap- 
proximations which guarantee gauge independence. These amount to com- 
pute the self-energies in (resummed) perturbation theory, satisfying the self- 
consistency condition in a perturbative sense. On the other hand, we shall 
not expand out the self-energies from eqs. (5) and (6) so that the resulting 
approximations for the S and M remain nonperturbative in g. This is possi- 
ble without encountering ultraviolet divergences because (in contrast to the 
original thermodynamic potential) the expressions (5) and (6) have the deci- 
sive advantage of being manifestly ultraviolet finite since the derivatives of the 
distribution functions vanish exponentially for both uj —>■ ±cxd. Moreover, any 
multiplicative renormalization with real Z leaves (5) and (6) unchanged. 



^ Vertices can be dressed in a self-consistent manner in the formalism worked out 
inRef. [17]. 



3 Contact with perturbation theory 



In Coulomb gauge, the gluon propagator consists of a transverse and a longi- 
tudinal piece, with LO self-energies at soft momenta {uj, k -C max(T, /i)) given 
by the so-called hard thermal (dense) loops [11,12] 
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is the (leading-order) Debye screening mass.[3]As is well known [18], the gluon 
propagator dressed by these self-energies has quasiparticle poles at UT^hik) 
with momentum- dependent effective masses and Landau damping cuts for 
|u;| < k. When k ^ m^, the pole corresponding to the collective longitudi- 
nal excitation has exponentially vanishing residue [19], whereas that of the 
transverse excitations tend to Jk"^ + m^, where the asymptotic mass is given 
by 



m^ 
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(9) 



The (massless) quark propagator at finite temperature or density is split into 
two separate branches of opposite ratio of chirality over helicity with prop- 
agators A± = [—u + fc ± S±]~^. In the HTL approximation, the respective 
self-energies read [18,12]: 



t±{u,k) 
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where {Cf = {N^ - 1)/2N) 



M' 
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^ For simplicity we sliall write out our formulae for only one value of the 
chemical potential; the generalization to several different chemical potentials is 
straightforward. 



The dressed fermion propagators have quasiparticle poles at uj±{k) and Lan- 



dau damping cuts. At large k and positive frequency lu^ -^ Jk"^ + M^ with 

asymptotic mass M^ = 2M^; the abnormal branch u!_ has exponentially van- 
ishing residue [19]. 

Although the various HTL self-energies constitute a leading-order result only 
under the condition of soft momenta and frequencies, the result for the asymp- 
totic masses applies also for hard momenta [20]. 

The LO interaction contribution oc g^ in our expressions for the entropy and 
the density arise from the domain of hard momenta. Because of the subtraction 
terms in (5) and (6), it turns out [14,16] that this contribution is given entirely 
in terms of the LO values for Ren7^(ct;^ = fc^) and IieT,±{uj = ±k), yielding 

9 Q ^6 -^2712 ^ ' 



{Ng = N"^ — 1) in agreement with known results [1].Q Notice the simplicity 
of these results in comparison with conventional two-loop calculations. Quite 
generally the LO interaction results for entropy and fermion density are given 
by 5(2) = -TI^Eb^Lb + ^EfM^^}, Ar(2) = -^J2^^pMlp where 
the sums run over all bosonic (B) and fermionic (F) degrees of freedom. Thus 
the LO interaction terms in S and A/" are entirely given by the quasiparticle 
dispersion relations at large momenta. In fact, these formulae for tS'-^) and A/'^^-' 
are identical to those corresponding to massive non-interacting particles (with 
constant masses m = moo and Moo, respectively) when expanded to leading 
order in m'^. 

When T ^ rho, the NLO contributions are of order g^ and in the conventional 
perturbation theory for the pressure arise from the sum of (electro-static) ring 
diagrams|f] which yields P^^^ = j^Ngifi^T oc g^. 

If the propagators and self-energies in (5) and (6) are those of the HTL ap- 
proximation, the order g^ contributions in the corresponding expressions for 
S and A/", denoted 5htl and A/rtl, turn out to give only the part of the full 
result that is obtained by differentiating P(^) at fixed rfiD [14,16], AT^tl = 
and iSjjTL = 12^^9'^i)- The remaining contributions of order g^ come from 
NLO corrections to the asymptotic quasiparticle mass and are calculable in 
HTL perturbation theory [16]. These corrections are, in contrast to the LO 
asymptotic masses, momentum-dependent, and increasing at hard momenta. 
In order to estimate their contribution in our numerical evaluations, we shall 



^ By contrast, the direct HTL resummation of the pressure overincludes the LO 

interaction contributions oc g^ [13,21]. 

"^ If r ^ g/j,, the ring diagrams give rise to contributions of order g"^ log{g). 



approximate them by averaged, effective NLO corrections to rrioo and Moo 
which are uniquely determined through eq. (12) as [16] 

5ml = -^g'NTmD, 5MI = -^g^CfTmo. (13) 



In these expressions, the hnear mo dependence and the Casimir factors are 
indeed those expected from the corresponding HTL-resummed one-loop dia- 
grams. 



4 Nonperturbative evaluation — recovering the pressure 



In our nonperturbative, numerical evaluation of (5) and (6) we consider two 
successive approximations. The first consists of using HTL-resummed propa- 
gators, which as we have discussed reproduces the LO interaction terms and 
part of the NLO one. As a second step we take into account the NLO cor- 
rections to the asymptotic mass in their averaged form (13) and we do so by 
^^L ~^ ""^oo ~ ""^oofl "" ^^^L/'"^oo]^^ ^"^^ similarly for the fermionic asymptotic 
mass Moo. This simple Fade resummation can be shown to give a surprisingly 
good approximation to the gap equation in exactly solvable scalar models 
[16]. However, because (13) refers to the asymptotic mass and because the 
plasma frequency and the Debye mass turn out to have rather different NLO 
corrections [22-24], we restrict this modification to hard momenta defined by 
k > A = v^SttTWdCa and vary c\ to test the stability of the final results. 

iShtl can be separated in two physically distinct contributions: one from quasi- 
particle poles, the other from Landau damping cuts of the various excitations. 
The gluonic ones, which only contribute to the entropy, have been given in 
eqs. (16) and (17) in Ref. [14]. The fermionic ones are quite analogous and 
read 5/htl = '^jhtl + '^/htl with 



S?W = 2iViV; l^^{T log(l + e-[-«')-'-l/n 

where only the explicit T dependences are differentiated and not those implicit 
in the dispersion laws uj+{k) and uj-.{k) of the fermionic quasiparticles. 

The contribution to the quark density is obtained by replacing d/dT in the 
above formula by d/d/j,. In the limit T — > 0, the resulting expression can be 



simplified to read (/i > 0) 

A^/i^TL ^^, = NNj I ^[^(^ - uj^ik)) - eiu^ik) - ^)]. (15) 



'' edk, 







The fermionic Landau- damping contribution to the entropy is 

S)l^^ = -mNj J ^^Oie - u;^){.rg[k - u; + S^^, k)] 

— Im S+(co', k)'Re[k — uj + S+(co', k)]^'^ 
+ aTg[k + uj + E4uj,k)]- lmJ:_{uj,k)Re[k + uj + J:^{uj,k)]^^} (16) 

and that to the quark density, A/'J'htl' is again obtained by replacing d/dT 
by d/d^. In the limit of zero temperature we simply have df /dfi —>■ 6{u — fi). 

The coupling constant in the HTL masses rh£, and M is determined from the 
two-loop renormalization group equation and as a central value for the (MS) 
renormalization scale fl we choose for the case of zero chemical potential the 
spacing of the Matsubara frequencies 2'irT, and for the case of zero temperature 
the diameter of the Fermi sphere, 2/i. The QCD scale Aj^ is, according to 
lattice data, rather close to Tc at /i = 0; for definiteness we adoptQ Tc/A^jg = 
1.14, a recent lattice result [25] for the case of pure glue QCD, which we take 
also for Nf j^ because lattice data so far showed little sensitivity of this 
parameter on A'^^. In order to have an estimate of the theoretical uncertainty 
of our approach we consider a variation of /2 around the central value by a 
factor of 2. 

Were it not for the trace anomaly [26] (T^) = S — 3P ^ the pressure would 
be simply given by P = {TS + /iA/')/4 (in the case of massless quarks). The 
correct relation is instead provided by an integration such as 

T 

P{T,fi) = Js{T',fi)dT' + P{T,,fi). (17) 



At yU = 0, the integration constant can be taken from the lattice; for definite- 
ness we take Ti = Tc and put P(Tc, 0) = 0. (Nonzero values as given e.g. by 
the lattice results do not modify the following results significantly for T > Tc.) 

In Fig. 1 our results for the entropy (a) and the pressure (b) of pure-glue QCD 



^ This value is about 10 % higher then the one used in Refs. [13,14]; the conse- 
quences for the following results are however rather small. 




Fig. 1. Comparison of our results (full lines LO; dash-dotted lines: NLO) for the 
entropy density (a) and the pressure (b) of a pure gluon plasma with the lattice 
results (grey bands) from Ref. [27]. See text for detailed explanations. 
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Fig. 2. Comparison of our results for the pressure for a quark- gluon plasma with 
two massless flavours with a continuum extrapolation of recent lattice results from 
Ref. [29]. 

are compared with the lattice results of Ref. [27][^. Full lines give the upper 
and lower bounds of the result in the HTL approximation under a variation 
of the renormalization scale as described above, the dash-dotted lines the 
corresponding one for the NLO approximation, which include a variation of 
ca in the range |. • -2. The lattice results are given by grey bands with the 
thickness giving a typical error. For comparison, the result of a direct HTL- 
resummation of the pressure as reported in Ref. [13], which fails to incorporate 
the correct LO term but reproduces the order ^f^-part, is given by the dotted 
line in Fig. lb (transcribed to our value of Tc/Aj^ = 1.14). 

Fig. 2 gives our results for the pressure in the presence of two massless quark 
flavours and compares with the (estimated) continuum extrapolation of recent 
lattice results [29]. The remarkably accurate agreement with the lattice result 
already at T > 2 . . . 3Tc is seen to persist also in the presence of fermions. 



° The more recent lattice results obtained with a RG-improved lattice action in 
Ref. [28] are consistent with the results of Ref. [27] within the given errors, but are 
systematically higher by some 3-4 %. 
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Fig. 3. Comparison of our results for the trace anomaly in a pure gluon plasma with 
the lattice result from Ref. [27]. 
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Fig. 4. Our result for the baryon density for Nf = 3 in comparison with order-g^ 
(dashed line) and order-(7*^ (dotted line) perturbation theory. 

In Fig. 3 the result for the trace anomaly is shown for the case of pure-glue 
QCD, in comparison with the lattice data from Ref. [27]. Here our results are 
very weakly dependent on /2 and also LO and NLO results almost coincide; in 
fact, at small T our results are dominated by the integration constant P{Tc, 0). 
Evidently, the details very close to Tc are not reproduced, but a substantial 
part of the lattice result is seen to be accounted for in our approach. 

Note that it is not possible to determine a particular -P(Ti) from the require- 
ment that both p{as) = P{T)/ Pq{T) and s{as) = S(T)/So(T) approach 1 in 
the limit of as — > 0, because the differential equation p{as) + ^P{as)p'{as) = 
s{as) with P{as) = —Poa'i — Pia^ — . . . has as homogeneous solutions 



p{a 



s/hoin. 



Ce 



-^[4/3(7'+0("«)] 



and this vanishes at a^ = together with all its derivatives! Different P{Tc, 0) 
thus correspond to exactly the same coefficients of a series expansion of p in 
as- Because it does not alter the entropy, the physical interpretation of this 
additional input is obviously that of fixing a (strictly nonperturbative) bag 
constant. 

At zero temperature we are lacking the information on the value of a critical 
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chemical potential fic and of P{^c)- However, we can still calculate A/", and 
the result of a numerical evaluation for iV = 3, Nf = 3 is displayed in Fig. 4. 
(NLO corrections are not needed now because the (/^-contribution is zero and 
therefore the HTL result is trivially correct at this order, too.) Also given 
are the perturbative results to order g^ and g'^ from Ref. [1,17] (translated to 
the MS-scheme). A simple quasiparticle model with constant mass m = M^o 
would result in values for A/" only slightly (~ 1%) above the perturbative 
gf^-result. Our results are still close to the latter, but somewhat higher. The 
perturbative order-g'^ results on the other hand are significantly lower with 
little overlap with the order-gf^ results. From the good agreement of our above 
results for the entropy with lattice results we expect that our approximations 
for the density are similarly stable for fi sufficiently above yUc- 



5 Outlook 



In our numerical evaluations, we have included NLO effects through the aver- 
aged values (13), which suffices to restore the correct coefficient of g"^ con- 
tributions. A more complete treatment would involve the exact result for 
SIIt{uj'^ = k"^) and 5T,±{ijj = ±k) in HTL-resummed perturbation theory. 
Their numerical evaluation and inclusion in S and A/" is work in progress. 

Another extension of the present results which is under way is the evaluation 
for general /x > and T > 0. There the Maxwell relations 

dS_ _dAf 



dT 



are fulfilled up to and including order g^ upon inclusion of the NLO contribu- 
tions. At and beyond order g'^, (18) constitute a nontrivial constraint on the 
renormalization-group flow of q;s(T, /i), determining P{T, fi) in terms of the ini- 
tial data at /i = 0, which can be taken from the lattice. This fixes the equation 
of state in the entire fi-T plane. A similar program has been carried out re- 
cently in Ref. [30] in simple quasiparticle models with momentum-independent 
transverse-gluon and quark masses. Within our approach, this can be extended 
to include both the effects of momentum- dependence of the thermal masses 
as well as Landau damping nonperturbatively, while maintaining equivalence 
with conventional perturbation theory up to and including order g^ in the 
thermodynamic potentials. 
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